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ABSTRACT 

Helicity and a effect driven by the nonaxisynimetric Tayler instability of toroidal 
magnetic fields in stellar radiation zones are computed. In the linear approximation 
a purely toroidal field always excites pairs of modes with identical growth rates but 
with opposite helicity so that the net helicity vanishes. If the magnetic background 
field has a helical structure by an extra (weak) poloidal component then one of the 
modes dominates producing a net kinetic helicity anticorrelated to the current helicity 
of the background field. 

The mean electromotive force is computed with the result that the a effect by the most 
rapidly growing mode has the same sign as the current helicity of the background field. 
The a effect is found as too small to drive an a^ dynamo but the excitation conditions 
for an aO dynamo can be fulfilled for weak poloidal fields. Moreover, if the dynamo 
produces its own a effect by the magnetic instability then problems with its sign do 
not arise. For all cases, however, the a effect shows an extremely strong concentration 
to the poles so that a possible aO dynamo might only work at the polar regions. Hence, 
the results of our linear theory lead to a new topological problem for the existence 
of large-scale dynamos in stellar radiation zones on the basis of the current-driven 
instability of toroidal fields. 



Key words: magnetohydrodynamics (MHD) 
- stars: interiors. 



instabilities - stars: magnetic fields 



1 INTRODUCTION 

Hydromagnetic dynamos can be understood as magnetic in- 
stabilities driven by a special flow pattern in fluid conduc- 
tors. There are, however, strong restrictions on the charac- 
teristics of such flows (see Dudley & James 1989 ) as well 
as on the geometry of the resulting magnetic fields (| Cowling 
I1933I ). The restrictions even exclude any dynamo activity for 
a number of flows. We mention as an example that differen- 
tial rotation alone can never maintain a dynamo (Elsasser 
1946). 

An open question is whether magnetic instabilities are 
able to excite a sufficiently complicated motion that together 
with a (given) b a ckgro und flow can generate magnetic flelds. 
iTout fc Pringlj ((1993) suggested that nonuniformly rotat- 
ing disks can produce a dynamo when magnetorotational 
(MRI) and magnetic buoyancy instabilities are active. Later 
on, n umerical simulations of Brandenburg ct al. ( 199a) and 
iHawlev. Gammie fc Balbuj (| 19961 ) have shown that MRI 
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alone may be suflicient for the accretion disk dynamo. It 
remains, however, to check (at least for the case of low 
magnetic Prandtl number) whether the MRI-dynamo has 
physical or numeric al origin (|Fromang fc Papaloizoul 120071 : 
iFromang erai]|2007l ). 

Another possibility was discussed bv lSpruitI l|2002l ) who 
suggested t hat different ial rotation and magnetic kink-type 
instability (|Tavleiil973h can jointly drive a dynamo in stellar 
radiation zones. The dynamo if real would be very important 
for the angular momentum transport in stars and their sec- 
ular evolution. It taps energy from differential rotation thus 
reducing the rotational shear. Radial displacements convert- 
ing toroidal magnetic fleld into poloidal fleld are necessary 
for the dynamo. The dynamo, therefore, unavoidably mixes 
chemical species in stellar interiors that may have observable 
consequences for stellar evolution. 

Such a dynamo, however, has not yet been demon- 
strated to exist. The doubts especially concern the kink- 
type instability that in contrast to MRI exists also without 
differential rotation. The Tayler instability develops in ex- 
pense of magnetic energy. Estimations of dynamo parame- 
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ters are thus necessary to assess the dynamo-effectiveness of 
this magnetic instability. 

The basic role in turbulent dynamos plays the ability of 
correlated magnetic (b) fluctuations and velocity (n) fluctu- 
ations to produce a mean electromotive force along the back- 
ground magnetic field B and also along the electric current 
J, i.e. 



{u X b) = aB — ^otjtJ- 



(1) 



We estimate the a effect by Tayler instability in the present 
paper. We do also find indications for the appearance of the 
turbulent diffusivity 77T in the calculations but we do not 
follow them here in detail. For purely toroidal fields we did 
not find indication for the existence of the term fix J which 
can appear in the expression ^ in form of a rotationally 
induced anisotropy of the diffusivity tensor. 

The fluctuating flelds for the most rapidly growing 
eigenmodes and the azimuthal averaging are applied in the 
LHS of Eq. IT} to estimate the a effect and its relation to 
the kinetic and magnetic helicity 'H'"" and H'^^'^'^. Our linear 
stability computations do not allow the evaluation of the a 
effect amplitude but its latitudinal proflle and its ratio to 
the product of rms values of u and b (i.e. the correlation 
coefficient) can be found. As the differential rotation is nec- 
essary for dynamo, we estimate also the influence of differ- 
ential rotation on Tayler instability. Next, a dynamo model 
with the parameters estimated for the magnetic instability 
is designed to find the global modes of the instability-driven 
dynamo. 



2 INSTABILITY AND HELICITY 

The model and the stability analysis of this paper are very 
close to that of iKitchatinov fc Riidiger! (12008) and will be 
discussed here only briefly. 

The basic component of the magnetic field inside a star 
is normally assumed to be the toroidal one. This toroidal 
field can be produced by differential rotation from even a 
small poloidal field. The background toroidal fleld of our 
model consists of two latitudinal belts of opposite polarities, 



r sin cos Oy/JIop f^Ae^ 



(2) 



(see Spruit 1999) with Q,a as the Alfven frequency of the 
toroidal field. Spherical coordinates are used with the axis 
of rotation as the polar axis and e^ as the azimuthal unit 
vector. The latitudinal profile of ((2)| peaks in mid-latitudes 
at 6* = 45° and 6* = 135°. 

The background fiow is simply 



U — r sin 9 Q e. 







(3) 



with n as the equatorial rotation rate. 

The f^A and Q are radius-dependent but this depen- 
dence is not of basic importance for the stability analysis. 
The reason is that the stratification of the radiative core is 
stable with positive 

9 dS 



N 



Cp dr' 



(4) 



where S is the entropy and Cp is the specific heat at constant 
pressure. 
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Figure 1. The ratio N/Qq in the radiation zones of the Sun (left) 
and of a 100 Myr old 3Mq star rotating with 10 days period 
(right). The model of the massive star is computed with the EZ 
code of lPaxtonI \2004 ) for X = 0.7, Z = 0.02. 



The buoyancy frequency TV is large compared to D, 
(N/Q ~ 400 in the upper radiative core of the Sun). Then 
the radial scale of unstable disturbances is short and the de- 
pendence of the disturbances on radius can be treated in a 
local approximation, i.e. in the form of exp(ifcr). The param- 
eter controlling the stratification infiuence on the instability 
is 

N 

Qkr' 



A = 



(5) 



and the most unstable d isturbances have A < 1 
(jKitchatinov fc Riidigeil 120081 ). This means that the radial 
scale of the disturbances, A = ir/k, is short compared to 
the radial scale SR of toroidal field or angular velocity vari- 
ations. In the solar tachocline where SI strongly varies in 
radius, the scale ratio is sma ller than unity, X/5R ~ 0.2 
(jKitchatinov fc Riidigen |2009| ). For such small scale ratio 
the radial derivatives in the linear stability equations are 
absorbed by the disturbances so that the local approxima- 
tion in the radial coordinate can be applied. Note that the 
unstable modes remain global in horizontal dimensions. 



2.1 Instability of toroidal fields 

The equation system of the linear stability analysis include 

^ + ([/ ■ V) u + (u ■ V) [/ + — (V (B ■ b) 
at flop 



- {B■\/)b-{b■\/)B) = - 



-VP 



P 



uAu 



(6) 



as the momentum equation for the velocity fiuctuations and 
the equations 



— =Vx(C7xb + uxB-?7Vxb), 



(7) 



for the induction of the magnetic fluctuations and for the 
disturbances (s) of the entropy 



ds 



CpX 



^ +U-Vs + u-VS^^:^^^AT'. (8) 

at T 

The equations ©-(O were reformulated in terms of scalar 
potentials for toroidal and poloidal parts of the magnetic 
and velocity flelds to reduce the number of equations, i.e. 

u — — curl((rVF)r) — curl curl((rl/)r) 
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-cui\{{rB)r) — curl curl((rA)r) 



( Chandrasekhar 
five eisenvalue 



1961). The 
equat i ons 



(9) 

resulting system of 

be found elsewhere 

(JKitchatinov fc Riidigen 120081 ) and it will here only be 
extended to the application of helical background fields. 
Again the eigenvalue problem is solved numerically. 

The equations include finite diffusion in opposition to 
the otherwise similar equations of Cally (2003) . The thermal 
diffusion is especially important because of its destabilizing 
effect. The Tayler instability requires radial displacements. 
It does not exist in 2D case of strictly horizontal motion. 
The radial displacements in radiation zones are opposed by 
buoyancy. The thermal diffusion smooths out the entropy 
disturbances to reduce the effect of stable stratification. This 
largely increases the growth rates for the instability of no t 
too strong, Q.a < ^, fields (jRiidiger fc Kitchatinovll201(J ). 
The diffusivities enter the normalized equation via parame- 
ters 












(10) 



In our computations we used the values e^ — lO"**, e,; = 4 x 
10~*, and Ely = 2 X 10"^° characteristic for the upper part of 
the solar radiation zone. As it must be the magnetic Prandtl 
number (5 ■ 10~^) exceeds the ordinary Prandtl number (2 ■ 
10-6). 

The stability problem allows two types of equatorial 
symmetry of unstable eigenmodes. We use the notations 
Sm (symmetric mode with azimuthal wave number m) and 
Am (antisymmetric mode) for these types of symmetry. Sm 
modes have vector field b which is mirror-symmetric about 
equatorial plane (symmetric hr,b^ and antisymmetric bg) 
and mirror-antisymmetric flow u (symmetric ue and anti- 
symmetric Ur,u^). Am modes have antisymmetric b and 
symmetric u. Instability can only be found for nonaxisym- 
metric disturbance s with m = ±1 in agreem e nt wi th the 
stability criteria of iGoossens. Biront fc Tavlen l|l98ll ) when 
applied to the toroidal field model of Eq. Q. 

The perturbations are considered as Fourier modes in 
time, in azimuth and radius in the form exp{i{kr+m(j}—ujt)). 
Only the highest-order terms in kr are used so that in ra- 
dial direction the theory is a local one (short-wave approx- 
imation). The wave number k enters the equations in the 
normalized form 



A 



N 
flokr 



(11) 



as a ratio of two large numbers. If only toroidal fields are 
considered with ^Ia/^ = 0.1 Kitchatinov & Riidiger (2008) 
found maximal growth rates of order 10"^ (normalized with 
the rotation rate) at radial scales of A ~ 0.1 for SI modes. 

As described in KR08, the eigenvalue u possesses a pos- 
itive imaginary part for an instability (the growth rate is 
7 = Sy(a;)). The equations have only be solved for the non- 
axisymmetric ('kink') modes with m, — ±1. One can easily 
show that the equation system possesses a symmetry with 
respect to the change of sign of the azimuthal wave number 
m. The equations are invariant under the transformation 



{m,uj,W,V,B,A,S)^ 

^ {~m,-Lj*,-W'',V,~B',A*,-S') 
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Figure 2. The normalized growth rate 7 = Qui/fl of S modes as 
function of the amplitude of the toroidal background field. Rigid 
rotation, m = ±1. Note the slowness of the modes for weak fields 
with Qa < ^ and the jump at Ha — H. The wave numbers 
are always optimized with respect to the maximum growth rates. 
Pm = 5- 10-3, Pr = 2- IQ-'^. 



where the asterisks mean the complex conjugate. If an eigen- 
mode exists for m = 1 with a certain growth rate 7 then the 
same 7 holds for the eigenmode with m, = — 1. We shall 
demonstrate the importance of this finding for the genera- 
tion of helicity and a effect. 

Figure [2] shows the growth rates for S±l modes in 
dependence of field strength. The results for A±l modes 
are very similar. For weak fields, 0.01 < f^A/fi < 1 the 
growth rates are closely reproduced by the parabolic law 
7 ~ Clf^A/^- For strong fields, Qa > fi, they are propor- 
tional to the field strength, 7 ~ SIa and do not depend on 
the rotation rate. 



2.2 Helicity production 

All averaging procedures in the present paper are realized 
by integration over the azimuth coordinate. Consider the 
formation of the kinetic helicity 



T-L '" = {u ■ curl u) , 



(13) 



(only the real parts of both factors) whose definition depends 
on the handedness of the used coordinate system. We shall 
prefer the righthand system. Let the expressions 



W 



(tor -I- iwije 



V ■ 



(vB. + \vi)e 



(14) 



represent the potential functions W and V . Then the real 
part of (I13p after the integration over the azimuth results to 



-K" 



sin^e 



mwi + sin 



de 



mvi + sin t 



^dwi 



+ 



+ mwji + sin 



-I dvK \ 



mv-n + sin 8 



dwR \ 



de J \ de J 



(15) 



(12) 



It is easy to show with the transformation rules ([12} that this 
expression for the modes with negative m has the opposite 
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sign as for the modes with positive m. Hence, 



■H'''"(m = -1) 



-V^''\m = 1) 



(16) 



(see Fig. [31 bottom). It means that for every unstable mode 
with finite helicity there is another unstable mode with the 
same growth and drift rates but with opposite helicity so 
that the resulting net helicity should vanish. Obviously, if 
all modes are excited the instability of a purely toroidal ax- 
isymmetric field can not produce finite values of the kinetic 
helicity. The same argument leads to the same conclusion 
for the current helicity 



H" 



— {b ■ curl b) . 



(17) 



There is also a more straightforward argument with the same 
result. The mode m — ~1 is identical to the mode m = 1 
but considered in a lefthand coordinate system. The sign of 
the helicity is equal in both lefthand systems and righthand 
systems. Hence, the mode m — —1 which gives the same 
helicity in the lefthand system as the mode m = 1 in the 
righthand system yields a negative helicity in the righthand 
system if the mode m = 1 yields a positive helicity in the 
righthand system. The net helicity in both the righthand 
system and the lefthand system, therefore, vanishes . Again 
the same is true for the current helicity (|17[) . 



3 INFLUENCE OF POLOIDAL FIELDS 

We find (see, however, Cally 2003) that from symmetry rea- 
sons unstable purely toroidal fields do not produce a net he- 
licity. An additional poloidal component of the background 
field, however, breaks the symmetry as then, for example, 
the two modes possess two different growth rates so that a 
certain sign of helicity will be preferred (Gellert, Riidiger & 
Hollerbach 2011). 

In the short-wave approximation only the radial com- 
ponent of the poloidal field is important. We write 



Br 



(18) 



with fj, = cos 9. Using the notation by KR08 with the oper- 
ator 

r=J^4sine4 + -4.^, (19) 



sin 6* 961 ""'"961 ' sin^ 6* 902 ' 
the equation for the azimuthal fiow reads 

uj (CW) = -i-^ (cw) + mCl (cw) 



mQ.A[CB) 



^W^|i((l-M')^)+-B|^((l-M^)f2A) 



+ {CV] 



d_ 

dfi 



((1 



lie 



d_ 

dfjb 



,2 do. 



CA] 



d_ 
9/i 



((I-m')^a) 



a^)-2(i 



dV 



Mr|f)-2(i-.^)^.) 



dA 
dfi 



^(r, trD\,n 2sdQA,p dB SOa.p . 



(20) 



The equation includes the poloidal background field via its 
last line. The parameter measuring the effect of the poloidal 
field is 
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Figure 3. Purely toroidal fields for rigid rotation. Top: The 
growth rate 7 for Qa = 0.3 vs. the radial scale A of the perturba- 
tions ylitl and Sitl. Bottom: The normalized small-scale helicity 
l|13|l . Solid line: m = 1, dashed line m = — 1. The net helicity van- 
ishes if both modes are simultaneously excited. Br = 0, A = 0.1. 



Oa.p 



iVnA,p _ iVB,.(/i) 



a n 



(21) 



Equation (|21|) shows that the characteristic strength of the 
field that can influence the Tayler instability (f^A.p ~ 1) is 
Qo/N times the toroidal field amplitude. This factor is of 
the order lO"'^' ■^'. 

For the latitudinal profile of the radial field component 
the simplest choice, i.e. Br(i-i) ~ /i, is used. Hence, both 
the background field components Br and B^ are antisym- 
metric with respect to the equator. The large-scale current 
helicity B ■ curlB, therefore, is also antisymmetric with re- 
spect to the equator. For positive amplitudes Ha and Qa.p 
it is positive at the northern hemisphere and negative at the 
southern hemisphere (it runs with cos'^ 0). We shall see that 
the pseudoscalar B ■ curlB alone determines the behavior of 
the pseudoscalars 'H'^'" and also of 'H'^"''. The basic rotation 
which via g ■ fl can also form a pseudoscalar does not play 
here an important role. 

The ratio /3 of the toroidal field amplitude and the radial 
field amplitude is 



B^ _ flA N 

Br fiA,p " ■ 



(22) 
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Figure 4. Helical background field. Latitudinal profile of the he- 
licity by the unstable modes Sitl (top) and Aibl (middle) for 
Qa = 0.3 and Qa.p = 0.1. Solid line: m = 1, dashed line: m = —1. 
The net helicity does no longer vanish if both modes are excited. 
It is ^ = 3N/n. Bottom: the S±l modes for U = 0. 



The equation for the meridional flow is 
uj{CV) 



-P (CS) - i^ (cv) 



A2 
- / - N / ON d(^iQ) dW odd 

^2,(1. (CB) +2(1-,^) ^1^^ +2m^^5 



-mClA (CA) - 2m 



djj, 



I — 2m (l~fi 



dfj, 
dClAdA 



+mn (cv) + 2m ^ V + 2m (l - /j,^) ^^^ 



l(nA.p{CA) + {l-fi') 



2,dQA,,dA ^90^^,(23^ 



dfi 9/i 



dfj. 



The equations for the magnetic fields components are 
uj (lb) = -i-^ (lb) + mC (QB) - ml (QaW) 



-m''—A~ — 
dn dfi 

2 OQa t r d 

Ofl Ofl 



2\2 dCl dA 
' dfj, dfj, 



2^2 dQAdV 
dfi dfi 



M 



^1{CiaA^w)hi-,, a^ Ofi 



^ m — T-^^V 



dfi 



and 



^{CA) = 



-i-^ (CA) + m(l (CA) - mClA {tv) 



The entropy equation is not infiuenced by the poloidal 
field. 

We computed the helicity "H"^'" of the critical modes 
always for the amplitude 17 a = 0.3 and CIa-p = 0.1 with and 
without basic rotation. The results are given in Fig. |4] 

The solid (dashed) lines give the helicity profiles for 
jn = 1 (m = — 1). Indeed, the helicity of the single modes 
is antisymmetric with respect to the equator (as also the 
current helicity of the background field). The plot at the 
bottom of Fig. |3]for the two modes SI is for 17 = 0. Note 
the extremely small differences to the top plot for the same 
modes under the infiuence of rotation. It is obviously the 
pseudoscalar B ■ curlB directing the formation of the helici- 
ties rather than the pseudoscalar g ■ fi formed by the global 
rotation. 

Figure [S] shows the normalized growth rates 7 = 7/t7 
vs the normalized wave length A for three different poloidal 
field amplitudes. From top to bottom: f7A,p = 0.3, f^A.p = 



0.1, f7A,, 



0.05. The peaks for the modes with m 



drift from A = 0.45 to A = 0.015. The corresponding peak 
values of the growth rates dramatically grow from 10~^ to 
10~^ demonstrating the strong stabilization of the Tayler 
instability for helical fields. An increase of the poloidal field 
amplitude by a factor of 6 (from bottom to top) leads to a 
reduction of the growth rate by three orders of magnitude. 
The main information of Fig. [S] is that the growth rates 
of the modes with m = 1 and tti = — 1 strongly differ. Ob- 
viously, they always produce helicity of opposite signs and 
with different growth rates. The growth rate for the SI mode 
with m = — 1 exceeds the growth rate of m = 1 by a factor of 
four. It is, however, much smaller than the rotation rate. The 
unstable Tayler modes are thus very slow, their correspond- 
ing growth times are much longer than the rotation periods. 
The helicity by the m, = — 1 mode is negative at the north 
pole opposite to the current helicity of the background field. 
The small-scale kinetic helicity and the large-scale current 
helicity are anticorrelated (cf. Gellert, Riidiger & HoUerbach 
2011, for a similar result in cylinder symmetry). 



4 THE ALPHA-EFFECT 

Linear stability computations do not provide the absolute 
value of a. Only its latitudinal profile and its relative mag- 
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Figure 5. The normalized growtii rates 7 for tiie modes Aibl 
and Sibl for Qa = 0.3 vs. the normalized wavelength A. From 
top to bottom: f2A,p = 0.3, ^A,-p = O-li f^A.p = 0.05. The modes 
with m = — 1 (solid line: S-mode, dashed line: A-niode) possess 
higher growth rates than the modes with m = 1 (dash-dotted: 
S-mode, dotted A- mode). From top to bottom the growth rates 
strongly increase while the characteristic length scale decreases. 
The product 'yX^ has always the value of about 10~®. 



nitude can be evaluated. We computed the normalized elec- 
tromotive force 



£ = 



N {ux b)4, 



(25) 



rmsi^rms 



which in opposition to the helicity (Fig.[4| is symmetric with 
respect to the equator. Urms is the rms velocity fluctuation 
after horizontal averaging (longitude and latitude), 6rms is 
the rms magnetic fluctuation. The factor N /Q, is introduced 
because the horizontal (velocity and magnetic) fluctuations 
are larger than the radial fluctuations by just this factor. 
With these normalizations the expression (|25[l is of order 
unity independent of the actual value of N/^. 
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Figure 6. The same as in Fig. [3] (bottom) but for the normalized 
electromotive force after Eq. I|25| l for the modes with m, = 1 (solid) 
and m = —1 (dashed). Ha.p = 0. 



Figure [6] gives the main results for purely toroidal fields 
with ([2]). The azimuthal component of the electromotive 
force i) vanishes at the equator, ii) the profile is symmetric 
with respect to the equator and iii) it is highly concentrated 
to the poles. The first finding is naturally for the a effect 
but it is not for the term f2 x J which could appear in the 
expression of the turbulence-induced electromotive force as 
a consequence of a rotationally induced anisotropy of the 
diffusivity tensor. We find, therefore, this effect not existing 
due to the Tayler instability of toroidal fields. For more com- 
plex field pattern its existence of cannot be excluded but it 
remained small in any case. 

For purely toroidal fields even the a effect does not exist 
as the modes with opposite sign of m (which have the same 
growth rates) do cancel each other not only with respect to 
their helicities but also with respect to the resulting EMF. 
In the nonlinear regime a spontaneous parity breaking may 
happen as it has been described by Chatterjee et al. (2010) 
and by Gellert et al. (2011). In this case, however, it might 
be impossible to predict the sign and the amplitude of the 
a effect. 

The concentration of helicity and EMF towards the 
poles reflects a basic property of the Tayler instability, i.e. 
that the instability pattern is more present in polar regions 
rather than in equatorial regions (Spruit 1999; Gaily 2003). 

Figure [S] shows the normalized growth rates 7 as a func- 
tion of the radial scale A under the influence of poloidal 
field components. The plots shows drastic differences of the 
growth rates between the modes of m = 1 (dotted lines) and 
m = —1 (solid lines). E.g., for C1a,p ~ 0.1 both modes with 
771 = —1 possess the maximum growth rates 7 with 5 • 10~^ 
at a wavelength of A = 0.15. The figure shows a very strong 
influence of the poloidal field amplitude on the growth rates. 
The growth rates are small for strong poloidal field but they 
are large for weak poloidal field. We are thus confronted 
with the dilemma that only background fields with finite 
current helicity originate fluctuations with a effect but the 
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0.05. The modes with the fastest growth produce an a effect 
which is always positive (negative) for positive (negative) 
current helicity at the northern (southern) hemisphere of 
the background field. The result is an a effect anticorrelated 
with the small-scale helicity and positively correlated with 
the large-scale pseudoscalar B ■ curlB. Exactly the same 
relations have been derived by nonlinear simulations of the 
kink-type instability for an incompressible fluid in a cylindric 
setup by Gellert et al. (20fl). We are thus encouraged to 
favor the results for the modes with the highest growth rates. 

For very weak poloidal field (f^A.p = 0.01, not shown) 
the Q effect is already so small (and its sign fluctuates) that 
the f^A.p = 0.01 seems to form the lower limit of the helicity 
production by poloidal fields. 

In the weak-field regime {CIa < 1) the instability excites 
stronger magnetic fluctuations rather than flow fluctuations. 
The ratio of fluctuating Alfven velocity Urms = fcrms/^/ioP 
to Mrms results as iirms/iirms — 470 for f^A = 0.01 aud 
I'rms/wrms — 46 for Qa = 0.1. Heucc, we find 



"Wrms 



4.7 

^A 



(27) 



for weak fields. For strong fields {Qa > 1) the fluctuations 
become close to equipartition, i.e. Umis — lirms. 




-50 50 

LATITUDE 

Figure 7. The normalized a effect parameter A for the modes 
with the highest growth rates (m = — 1, identic for A mode and 
S mode). From top to bottom: f^A^p = 0.3, f^A^p = 0.1, CIa^p = 
0.05. Note i) the extreme concentration of A to the poles and 
ii) the A amplitudes strictly running with f^A.p- For the given 
magnetic geometry the a effect is always positive (negative) at 
the northern (southern) hemisphere. Qa = 0.3. 



corresponding poloidal field components suppress the Tayler 
instability (see Rfidiger, Schultz & Elstner 2011). 
For the further discussion the quantity 



A 



£ 



sin o cos f 



aBo N 



(26) 



.'rms'.'rins 



is introduced which is antisymmetric with respect to the 
equator as expected for the (normalized) a effect. The 
sin cos 6 in the denominator eliminates the latitudinal pro- 
file of the toroidal field and Bo — ^/JmiP rflA is its amplitude 
(seeEq. (0). 

Also the resulting ^-profiles are highly concentrated to 
the poles (Fig. [7}. The plots show the profiles for the weak- 
field cases Qa ~ 0.3 with CIa,p = 0.3, 17a, p = 0.1 and Ua.p = 



5 DYNAMO THEORY 

On the basis of the obtained informations about the a ef- 
fect (amplitude and latitudinal profile) we have to probe the 
possible existence of a dynamo mechanism. 



5.1 Eddy difFusivity 

We start to estimate the typical velocity perturbation and 
the eddy magnetic-diffusivity. We write u ~ X/t and use the 
growth time Tgr — I/7 as the timescale r. Hence, 



7A: 



-7A — . 



(28) 



The first step holds for purely toroidal fields (Oa.p ~ 0). 
From Fig. [3] (top) one finds 7A = 10~* so that with solar 
values M ~ 1 mm/s results. The estimate 77T — uX leads to 
77T — 10 cm /s as the order of magnitude of the eddy diffu- 
sivity. The value increases for fast rotating giants by one or 
two orders of magnitudes. It perfectly fits the diffusion coef- 
ficients for chemicals which are needed to explain the weak 
lithium depletion of solar type stars (Barnes, Charbonneau 
& MacGregor 1999). 

The same estimate rjr — 'yX? yields 



77T 



2- ?2 ^ o 2 



(29) 



The three examples for helical background fields given in 



Fig. [5] lead to a common value of 7 A 



10" 



Hence, r/x 



3 • 10^ cm^/s results for the upper radiation zone of the Sun. 
This value is very close to the above estimation for purely 
toroidal fields. It is increased by more than three orders of 
magnitudes if more massive stars are considered. Figure [T] 
shows the details of a model with a mass of 3M0 and 10 
days rotation period. 
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5.2 Dynamo conditions 

The a effect enters tfie dynamo theory via the dimensionless 
dynamo number 

C, = — . (30) 

Equation (|26|) together with the heuristic relation u^/r/x — 7 
directly provides 






(31) 



Hence, the product 7A determines the effectivity of the a 
effect. With 7 yl ~ 0.03 a value of Cq ~ lO"'^"^^ resuhs for 
the solar model. Note that the buoyancy frequency linearly 
enters the equation. Also for the ?>Mq star of Fig. [T] the 
Ca — 0.1 remains small. The operation of a (stationary) 
o? dynamo is thus excluded for radiation stellar zones. The 
existence of an (oscillating) afJ dynamo, however, remains 
possible. All oTl dynamos can work with very small a effect 
if only Cn 



Co 



!^0 T^ 



(32) 



is high enough and this is always possible for sufficiently 
small eddy diffusivity r/x. The only consequences of very 
small a effect are i) that the ratio of toroidal and poloidal 
magnetic field components becomes very large and ii) also 
the growth time of the dynamo instability becomes large. 
The growth time of the weakly supercritical aQ dynamo 
for rjT — lO^cm^/s is extremely long (order of Gyrs for the 
Sun). With (say) 3% differential rotation this value of the 
eddy diffusivity leads for the Sun to Cn of order 10^. 

There is another possibility to proceed. For all aQ dy- 
namos the ratio /3 of the toroidal and the radial field ampli- 
tudes is 



(33) 



The scaling parameter e is about 0.05 (see Table [T}. Hence, 
Cn ~ ff'Caj^ ■ The excitation condition for such dynamos 
can be written as 




(34) 



CaCn ~ Dc 



where Dcrit runs inversely with the shear. With (|33p 
Cq, > Ve^ -Dcrit//? results as excitation condition. One finds 



7A > -V-Dcrit^A.p^A 

5 



(35) 



for fixed shear. Note that the condition (|35p does no longer 
contain the stellar parameters like Q,/N and/or the radius. 
The following findings are thus valid for all stellar radiation 
zones. 

Our first example is formed by f^A = 0.3 and nA,p = 
0.3. The condition (PSj) then reads 7A > 0.06e\/D'crit. As we 
shall see below, a typical value for the dynamo number is 



: 30 so that for dynamo excitation 



^A > 2nA,p. 



(36) 



This relation must be read as an equation for the supercrit- 
ical value of r^A.p. The excitation is thus formally more easy 
for smaller poloidal field. After the numerical results in the 
Figs. [5] (top) and[7](top) one finds ^A ~ 10"^ for Cl a,p = 0.3 



which is by far too small. For f^A.p = 0.1 it is ^A ~ 0.03 
which is also not supercritical. The comparison of the Figs. 
[S] and [7] demonstrates that for large poloidal fields the A 
grows but in the same time the growth rate drastically sinks 
which indicates the stabilizing action of the poloidal fields. 
An af2 dynamo in radiation zones cannot work, therefore, 
with too strong poloidal fields. 

Note that with 'yA ~ 0.2 for Q,a,p ~ 0.05 the condition 
H36[) can be fulfilled due to the increase of the growth rate. 
For smaller fJA.p the ability of the poloidal field to create 
a coherent a effect sinks. For CIa,p ~ 0.01 the function A 
has already two signs at either hemisphere. However, as the 
condition (|36|) can be fulfilled the existence of an aSl dynamo 
cannot be excluded by the numerical results. We must thus 
solve the dynamo equations in order to probe the existence 
of such dynamos which work with very small values of the 
Q effect. 

The a effect plotted in Fig. [7| shows another compli- 
cating characteristics, i.e. its concentration at the poles. In 
the next Section nonlinear ail dynamo models are thus con- 
structed with a weak solar-type latitudinal differential rota- 
tion and an a effect which is concentrated at the poles. It is 
known from the theory of aQ dynamos operating with cos 9 
profiles of the a effect that they produce too polar butter- 
fly diagrams. We must thus expect that the existence of aQ 
dynamos to produce mid-latitude belts like in Eq. ^ with 
pole-concentrated a effect is unlikely. 



5.3 Dynamo models 

The influence of the polar concentration of magnetic- 
induced a effect on a possible ail dynamo can easily be 
probed. The stability of an axisymmetric toroidal back- 
ground field has been considered so that only the numer- 
ical values of an axisymmetric a are known. This scenario 
is only consistent if the possible ail dynamo produces ax- 
isymmetric toroidal field belts at the latitude of the original 
background field. We know that this is true for an a ef- 
fect (plus solar-type differential rotation) with the standard 
cos ^-profile. We have to check, therefore, whether this result 
remains true for the a effect concentrated close to the poles. 
To this end a simple model is constructed. In a spheri- 
cal shell between the normalized radii 0.6 and 1 the rotation 
frequency is f2 = Cn(l — 0.03cos^&). The small value of 
the relative shear allows to consider the fluid as hydrody- 
namically stable (see Watson 1981). Possible production of 
hydrodynamic-induced a effect in radiative zones such as 
that by Dikpati & Gilman (2001) is therefore excluded. The 
shear peaks at 45° which is the same latitude as that of the 
peak of the toroidal fleld. To model the polar concentration 
of the magnetic-induced a effect the expression 



a = Ca cos 



(37) 



is used with the free parameter £ fixing the latitudinal pro- 
file of the a effect. The polar concentration presented in 
Fig. [7] leads to rather high values of £. A perfect-conductor 
boundary condition is used at the inner radius. The outer 
computing domain is extended to the radius 1.2. Outside 
the stellar surface the diffusivity is increased by a factor of 
10. For the outermost boundary a pseudovacuum condition 
is used. The method of a global quenching of the a effect is 
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used in order to find tlie cliaracteristic eigenvalue of marginal 
dynamo instability (Elstner, Meinel & Riidiger 1989). 

Table [l] presents the numerical results for a reference 
dynamo model with 3% difTerential rotation and for Cn — 
257, 000. Only the numbers are given for marginal dynamo 
instability. If the real Ca differs from the reference Ca by a 
factor 1/7 then the resulting P also differs by the factor 7, 
so that 



/3Cc 



P''''C'f 



30. 



(38) 



This condition is automatically fulfilled for all dynamo mod- 
els which fulfill the excitation condition (|35p independent on 
the particular stellar model. All the dynamos listed in Ta- 
ble [T] provide positive (negative) values of the large-scale 
current helicity B ■ curlB at the northern (southern) hemi- 
sphere. As we have assumed the same constellation for the 
above background field producing the a effect the theory is 
thus consistent under this aspect. 

We have also to probe whether the star rotates fast 
enough to produce sufficiently large Cn for the required dy- 
namo numbers of the reference values of Table [T] It is 



Cq = 



Dc 



/D~r ^A N 



Lyct 



^A,p 



Q' 



(39) 



Note at first that the excitation is much more easy for mas- 
sive stars as their N/il is smaller (Fig. [T|). Inserting the char- 
acteristic numbers of our models the result is Cn ^ lO'' N/Q 
which is fulfilled by both the considered stellar models by 
their estimated values Cn ~ 10^'' . Either the Sun as the con- 
sidered hot stars with rotation periods of a couple of days 
rotate fast enough the excite an aQ dynamo with a very 
weak shear (3%). 

To discuss the infiuence of the quantity £ we present 
models with i — 1, £ = 5 and £ — 15. The latter value de- 
scribes the strongest concentration of the a effect to the 
poles. For £ = 1 (i.e. a (x cos 6) the dynamo-generated 
toroidal field peaks exactly at 6 — 45° so that the presented 
stability analysis which bases on a toroidal field which also 
peaks at 9 = 45° (see Eq. ^) would be consistent. This case 
together with the excitation condition p6|) which is fulfilled 
by the magnetic amplitudes CIa = 0.3 and f^A.p ~ 0.05 
would strongly suggest the existence of an aQ dynamo in 
magnetized radiation zones of hot stars. 

However, the models listed in Table [T] provide the re- 
sult that the induced toroidal field belts become more and 
more concentrated to the poles for increasing £. This be- 
havior is insofar not trivial as the generation of the toroidal 
field by the differential rotation does not depend on £. Such 
polar belts resulting for £ > 1 can never reproduce the mid- 
latitudinal belts of Eq. @ on which the a effect bases. The 
first example of Fig. [8] (top) is the standard ail dynamo 
with the latitudinal profile a oc cos 9. It produces axisym- 
metric toroidal magnetic belts of dipolar symmetry at the 
same latitude where the shear dQ,/d9 has a maximum. For 
growing £ the belt position drifts more and more polewards 
(the field position becomes 45° > 9 > 135°). It is thus not 
possible to maintain toroidal field belts in mid-latitudes if 
the a effect mainly exists in the polar region {£ > 1). 








0.0 0.^ O.'l 0.6 0.» 1.0 1.2 




0.0 0.2 0.1 O.ti O.H 1.0 





0.0 o.a 0.4 0.6 o.[i 1.0 i.a 



Figure 8. The field geometry of of2 dynamo models with 1=1 
(top), £ = 5 (middle), I = lb (bottom). Left: radial field at the 
surface. Right: toroidal field pattern. The magnetic field belts are 
drifting radially during the cycle. 
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Table 1. Reference dynamo models with 3% differential rotation 
and for Cq = 257, 000. The last column gives the polar distance 
of the maximum position of the northern toroidal field belt. 





^crit 


Cc 


/3 


e 






i 


eV-Dcrit 


e 


1 


118,000 


0.46 


56 


0.07 


26 


45° 


5 


347,000 


1.35 


23 


0.05 


31 


38° 


15 


976,000 


3.8 


6.3 


0.02 


24 


23° 



6 CONCLUSIONS 

It is demonstrated with a linear theory that the Tayler in- 
stability of a toroidal magnetic field in a density-stratified 
radiation zone of hot stars does not produce helicity and/or 
a effect. If, however, a weak poloidal field component is 
added forming a large-scale current helicity B ■ curlB then 
the instability leads to small-scale helicity, current helicity 
and also to a effect. If B ■ curlB is positive in the north- 
ern hemisphere then the helicity {u ■ curlu) is negative at 
the northern hemisphere if only the modes are considered 
which grow fastest. Hence, the small-scale kinetic helicity 
and the large-scale current helicity of the background field 
are anticorrelated. 

As it is often the case, the corresponding a effect is anti- 
correlated with the kinetic helicity and, therefore, positively 
correlated with the large-scale current helicity B ■ curlB. 
The calculations lead to two basic properties of this a ef- 
fect. It is i) small, i.e. the normalized value Ca is only of 
order 10~""^' and ii) concentrated at the poles. The first 
property excludes the existence of a^ dynamos in radiative 
zones and the second property makes the existence of afl 
dynamos (with weak differential rotation in mid-latitudes) 
unlikely. As we have shown only a dynamo with a ex cos^ 
produces the toroidal fields in mid-latitudes. For £ > 1, how- 
ever, the belts are more and more shifted into the polar re- 
gion. Such fields cannot close the loop of field amplification 
by reproducing the original axisymmetric toroidal field. 

Or, with other words, the rotation law with the consid- 
ered profile (n (X cos^ 9) mainly induces toroidal fields at 
mid-latitudes where for high values of £ almost no a effect 
exists. Hence, an aQ dynamo could only work for very fast 
rotation. By very fast rotation the instability is suppressed 
(see Fig. [2]). This topological problem seems to be the key 
problem with the magnetic-driven dynamo rather than the 
excitation conditions for aQ dynamos. It is, however, not yet 
clear whether this argumentation also holds with the same 
power with other than the used rotation laws and/or in fully 
nonlinear simulations. All the presented aSl dynamo mod- 
els working with a weak solar-type differential rotation are 
reproducing the assumed sign of the large-scale current he- 
licity B-curlB. If thus the dynamo produces its own a effect 
by the magnetic instability then the signs will be consistent. 

A basic deficit of the presented theory is the fact that 
both helicity and a effect have been computed in rigidly ro- 
tating stars while for an aQ dynamo the rotation must be 
nonrigid. It might be the case that growth rates and mag- 
netic patterns of the Tayler instability are strongly modified 
by even weak differential rotation. The discussion of this 
new subject, however, is not the scope of the present pa- 



per. We have thus considered only cases with a rather weak 
differential rotation (3 %). 

Another deficit is formed by the order-of-magnitude es- 
timation of the eddy diffusivity. Using the characteristic 
scales of the modes with the highest growth rates the typ- 
ical velocity is of order 1 mm/s and the typical diffusivity 
value is about 10^ cm^/s (both for solar values). Estimations 
of the radial mixing produced by the same instability pro- 
vide an important test of the theory. The observed content 
of fight elements at the Sun impose restrictions on radial 
mixing in the upper radiative core (Barnes et al. 1999, and 
references therein). The chemical mixing in the deep stel- 
lar interior can also be compatible with observations only if 
its ch aracteristic time is longer tha n the evolutionary time 
scale (jKippenhahn fc Weigertlll994l ). Only very slightly su- 
pcrcritica l kink-type instability can satisfy this restriction 
(Kitchati nov fc Riid igcr 2008). 

All the material values in the paper are solar values. 
The formulation of the dynamo theory is such that the ma- 
terial parameters of the stellar interior do not appear (see 
Eq. (|35|l l. The conclusions about the dynamo activity for hot 
stars do thus not depend on the mass of the star. This is not 
true, however, for the above mentioned order-of-magnitude 
estimations of the characteristic values of velocity and diffu- 
sivity which are running with fl/N and (fl/N)'^, resp. From 
Fig. [T] one finds that the given solar values increase by fac- 
tors of 30 or 900, resp., for stars with three solar masses. 

Only nonlinear simulations can demonstrate whether 
the system indeed prefers the modes with the highest growth 
rates. Only if not then the radiative-zone dynamo has a 
chance to work in the solar interior but if it works then 
it should be only marginally supercritical. 
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